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A possible sequential star formation in the giant molecular 
cloud G174+2.5 
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ABSTRACT 

We investigate the nature of 14 embedded clusters (ECs) related to a group of four 
H II regions Sli2-235, Sli2-233, Sh2-232, and Sh2-231 in the giant molecular cloud 
G174 + 2.5. Projected towards the Galactic anticentre, these objects are a possible 
example of the collect and collapse scenario. We derive astrophysical parameters {age, 
reddening, distance, core and cluster radii) for the ECs and investigate the relationship 
among their parameters. Parameters are derived with field decontaminated 2MASS 
colour-magnitude diagrams (CMDs) and stellar radial density profiles (RDPs). The 
CMDs of these young clusters are characterised by a poorly-populated main sequence 
and a significant number of pre-main sequence stars, affected by differential reddening. 
The ECs are KKC 11, FSR 784, Sh2-235 E2, Sh2-235 Cluster, Sh2-233SE Cluster, 
BDSB 73, Sh2-235B Cluster, BDSB 72, BDSB 71, Sh2-232 IR, PCS 2, and the newly 
found clusters CBB 1 and CBB 2. We were able to derive fundamental parameters 
for all ECs in the sample. Structural parameters are derived for FSR 784, Sh2-235 
Cluster and Sh2-235E2. 

Key words: (Galaxy.) open clusters and associations:general; Galaxy: stellar content, 
collect and collapse, and sequential star formation; Galaxy: structure 



1 INTRODUCTION 

Most stars form after the gravitational collapse of massive 
and dense gas clumps inside giant molecular clouds (GMCs), 
with collapsing clumps forming embedded clusters (ECs). 
However, supernova explosions, H II region expansion due to 
massive-stars, UV radiation, and stellar winds, may disrupt 



GMCs completely, on a ti mescale of a few 1 0^ vr ([Elmegreen 



l200d: iBonnell et all [200^ ). In this context. iHartmann etal 
( 200ll ) point out that stars older than ~ 5 Myr are not 



found associated with molecular gas, and lAllen et al.l (|2007l ') 
sugge st that the primordia l gas of ECs disperses in 3 — 5 
Myr jLeisawitz et al.l [19891 : iProszkow fc Adain^ 120091 '). For 
iLada fc Ladal (|2003h . the duration of the embedded phase 
is 2 — 3 Myr. On the other hand, t he parent GMCs are 
very disruptive environments for ECs. iLada fc Ladal (|2003l ) 
estimate that only 4 - 7% of them survive for more than 40 
Myr {infant mortality), but the appearance of bound and 
unbound clusters are indistinguishable for clusters younger 
than 10 Myr. In other words, stars are born in star clusters 
embedded in GMCs but, as a consequence of the disruptive 
mechanisms, most end up as part of the field-star population 
after cluster disruption. 

ECs can be partially or fully immersed in embryonic 
molecular clouds and H II regions. The youngster are 
located in gas clumps and the most evolved are often 



linked with H II regions or other nebulae (|Leisawitz et al.l 
|l989). Often, these H II regions, excited by fast winds 
from massive OB stars, expand into the molecular cloud 
triggering sequential star form ation. In this sen s e, the 
collect and collapse model (|Elmegreen fc Ladal Il977l : 
IWhitworth et al"]|l994l ) suggests that H II regions expand 
accumulating material between the ionisation and the 
shock fronts. This material, as a consequence of the 
shocks, becomes unstable and fragments into several 
cores, triggering star formation in multiple protocluster 
regions. Another possible process is the "radiation-driven 
implosion" model, in which the expanding H II region 
compresses the existing molecular clumps, the density 
increase exceeding t he cri tical mass and the clump collapses 
(Lcfloch fc La zarcfj Il99j ). In both scenarios the massive 
stars t rigger a second gen er ation of cluster formation (see 
also, iFukuda fc Hanawa' '2000'; 'Hos okawa fc Inutsukal 
20051: [Dcharvcng, Zavagno fc Caplm] l2005l : 



Dale. Bonnell fc Whitworthi 120071 ). In any case, GMC 



observations indicate that there often occurs a multiple 
cluster formation as a result of winds of OB stars or 
expanding H II regions lAUen et al.ll2007l ). 

Star cluster formation is the preferential mode of star 
formation, and ECs may be responsible for 70%-90% of all 
stars formed in GMCs. However, the rapid expulsion of the 
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Table 1. Cross-identifications of the embedded clusters in G174 + 2.5. 
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Figure 1. XDSS R image (30' X 30') of the Sh2-235 HII region. 
The image also shows part of the Sh2-232 HII region (top left). 
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primordial gas by winds of OB stars and supernova explo- 
sions disrupt most clusters very early. As a consequence of 
rapid gas expulsion, the stellar orbits cannot adjust to the 
new potential, and probably give rise to an unbound associ- 
ation. N-body simulations show that in this phase (10 - 30 
Myr) the cluster expands in all scales reaching for virializa- 
tion. The fate of a cluster is determined by the star formation 
efficiency (SFE). If the gas is removed slowly the cluster will 
remain bound as long as the SFE is higher than 30%, but if 
the gas is removed rapidly, the SFE needs to be higher than 
50%. As a consequence of infant mortality, the number of 
optically detected clusters is significantly srnaller than that 
of ECs (iLada et al.lll984l: IVerschuerenI Il990l : iLada fc LadH 
l2003l : lGoodwin fc Bastianll2006l ). 

Young clusters in expanding H 11 regions are charac- 
terised by a significant population of low-mass pre-main se- 
quence (PMS) s tars together with variable fractions of mas - 
sive OB stars JSaurin et al.l |201GI : iBonatto fc Bical I2OIII ). 
iDobashi et al.l l|200ll ) show that protostars associated with H 
II regions are more luminous than those in molecular clouds 



away from them, which indicates that H 11 regions favour 
massive star s or cluster form a tion i n neighbouring molec- 
ular clouds. iDeharveng et al" I (|2003') suggest that clusters 
formed from these processes are still deeply embedded in H 
II regions, and are not dispersed. 

The observed colour magnitude diagram (CMD) of 
an EC in general presents a relatively vertical, poorly- 
populated main sequence (MS) and a well developed 
PMS with a large p opulation of faint and red stars 
(|Bonatto fc Bical l20od ). This feature is evident when the 
raw photometry and decontaminated CMDs are compared 
to field star extractions in the present analysis. The low 

mass stars rema in in the PMS phase for ~ 30 Myr 

(jGouhermis et al.ll2008l ). 

In this work we investigate properties of the star clus- 
ters embedded in the GMC G174-I-2.5, related to a group of 
H 11 regions located in the Perseus arm towards the Galac- 
tic anticentre. A star forming complex with numerous ECs 
minimises uncertainties owing to the distance in common. 
This paper is organised as follows. In Sect. [2] we provide 
general data on the target clusters. In Sect. [3] we present 
the 2MASS photometry and introduce tools employed in 
the CMD analyses, especially the field star decontamination 
algorithm. Sect. |4] is dedicated to the cluster structure. In 
Sect.[5]we estimate cluster masses. In Sect. |6]we investigate 
the relationship among the derived parameters. Finally, in 
Sect.[7]we present the concluding remarks. 



2 PREVIOUS STUDIES ON THE G174+2.5 
STELLAR CONTENT 

The H II regions Sh2-235, Sh2-233, Sh2-232, and Sh2-231 
have been subject to several studies in different wavelength 
ranges, but astrophysical parameters for most ECs located 
in these regions have not been derived so far. The ECs se- 
lected for the present analysis are given in Table (TJ where 
we adopt a chronological criterion for the literature identifi- 
cations. We show in Fig.[T]an XDSS 30' x 30' R image of the 
Sh2-235 H 11 region. In Fig. [2]we show schematic charts for 
the nebulosities and clusters, and in Fig. |3] 2MASS images 
of the ECs in the Ks band. All objects are quite obscured 
by interstellar dust and appear to be affected by differential 
reddening. Being projected close to each other, this sample 
of ECs may be an example of H 11 regions with a sequential 
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Figure 2. Top left panel: schematic map of the nebulae (open circles). Filled circles are ECs, asterisks are O stars, and plus signs are B 
stars. Top right panel: zoom of the region that contains the HII region Sh2-235. Bottom panel: zoom of the gray region in the left top 
panel that contains the nebulae Sh2-231, Sh2-233, RN049 and Sh2-233SE. 



star formation and the development of a collect and col- 
lapse scenario. Furthermore, the detection of far-IR sources, 
molecular outflows, and H2O masers in previous works in- 
dicate on-going star formation. The distance to the Sun es- 
tima ted for these H II regions is in the range (Iq = 1.0-2.3 
kpc l|Georgelinlll975l '). Most works use = 1.8 kpc for this 
group of nebulosities. 



Within uncertainties, Sh2-231, Sh2-232, Sh2-233 
and Sh2-235 have com parable CO radial velocities 
ijBlitz. Fich fc Stark|[l98^ ). This indicates that we are deal- 
ing with a large Hll-molecular complex with components 
located essentially at the same distance from the Sun. 



2.1 Sh2-235 and surroundings 

Sh2-235 is the most prominent H II region in this group. It 
is a diffuse optical H II region ex cited by a sta r of spectral 
type 09.5 V (BD+35°1201). Alle n et al.1 (|2005D identify two 
clusters associated with the Sh2- 235 H II regi o n, Sh 2-235 
Cluster and KKC 11. iXumar. Keto fc ClerkinI (I2OO6I) add 
to the se objects the cluster Sh2-235 East2. Ixirsanova et al.l 
(|2008l ) conclude that these objects are still embedded in 
dense clumps of the parental molecular cloud 6*174 + 2.5. 
They also argue that the Sh2-235 Cluster and Sh2-235 E2 
probably started the primordial gas expulsion, but KKC 11 
is less evolved. 

South-west of the Sh2-235 are four small nebulae named 
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Figure 3. First row panels: 2MASS 3' X 3' Ks images centred on FSR 784, KKC 11 and Sh2-235E2. Second row: the same for the 
cluster pairs Sh2-235 Cluster and CBB 2, G173 and CBB 1, PCS 2 and Sh2-233SE Cluster. Third row: Sh2-235B Cluster, BDSB 71, 
and BDSB 72. Fourth row: same as row three for BDSB 73 and Sh2-232 IR Cluster. The small circle is automatically generated by the 
2MASS image tool, surrounding the input coordinates. The large circle in the second row surrounds the main body of each cluster. The 
last panel in the fourth row shows the large background area (indicated with outer and inner brown circles) used in the statistical field 
decontamination of Sh2-232 IR Cluster. Note that it is huge with respect to the cluster area, for statistical purposes. The symbols are 
the same as in Fig. |2] 
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Figure 4. 2MASS CMDs extracted from the central regions of 
the cluster pair Sh2-235 Cluster and CBB 2. Top panels: ob- 
served CMDs J X (J — Ks). Middle panels: equal area comparison 
field. Bottom panels: field star decontaminated CMDs fitted with 
the Padova isochrone for MS stars and Siess for PMS stars. The 
colour-magnitude filters used to isolate cluster MS and PMS stars 
are shown as shaded regions. We also present the reddening vec- 
tor for Ay = to 5. BD-|-35°1201 is shown as a circle around the 
star in the top and bottom-left panels. 



Figure 5. 2MASS CMDs extracted from the central regions of 
FSR 784, KKC 11 and Sh2-235 E2, respectively. Top panels: ob- 
served CMDs J X (J — Ks). Middle panels: equal area comparison 
field. Bottom panels: field star deconta minated CMDs fitted with 
the Pa dova isochrone for MS stars and lSiess. Dufour &: Forestinil 
ll2000l 'l for PMS stars. The colour-magnitude filters used to isolate 
cluster MS and PMS stars are shown as shaded regions. We also 
present the reddening vector for Ay = to 5. 



GGD 6 (also known as RNO 52), Sh2-235 A (GM1-G6), 
Sh2-235 B (BPS 46 or GM1-G5) and Sh2-235 C (GGD 5 
or BFS 47). Sh2-235 A and Sh2-235 B are located about 
10' (5.6 p c) south of Sh2-235. and separated by 40" 
(0.35 pc) (jBlitz. Fich fc StarklflQsl ). Sh2-235 A is a com- 
pact H II r egion of » 2 " (0.1 7 pc) in diameter (jFelli et al.l 
liggTl ). For iBolev et al.l (|2009l ') the exciting star of Sh2-235 
B is an early- type Herbig Be star of spectral type BIV. 
iHodappI ^9ML found a cluster i n Sh2 - 235 B, which was 
confirmed bv'Wang. Testi fc Fellil (|l997l '). lAllen et all (|2005l ) 
and Kirsanova ct al. (2008). Sh2-235 C is a small H II reg ion 
located about 3.5' south of Sh2-235 B (jFeUi et al.ll2004i l. It 
coincides with an optical nebula and is excited by a BO. 5 
star, presenti ng a Herb i g-Har o object and a par tial shell 
morphology (|Felli et al.1 |l993). iBica et al.1 (120031 ) identify 
other three ECs in these nebulosities, BDSB 71, 72 and 73. 

Based on the spatia l distribution of c luster s and kine- 
matics of molecular gas iKirsanova et al.l ()2008l l point out, 
that the expansion of Sh2-235 would be responsible for clus- 
ter formation in that area. They suggest that sequential star 



formation is triggered by a combination between compres- 
sion of pre-existing dense clumps by the shock wave, and 
the collect and collapse scenario. However, the clusters in 
Sh2-235 A, B and C appear to be embedded in primordial 
gas, and star formation cannot be t riggered by the ex pan- 
sion of the Sh2 -235 ionisation front llLafon et al.lll983 ). On 
the other hand, iTokunaga fc Thompsonl "i 19791 ) suggest that 
the linear disposition of these star forming regions may be 
the result of collect and collapse scenario. 



2.2 Sh2-231, Sh2-232 and Sh2-233 

Sh2-232 is an extended H I I region with 40' diameter 
(Fig. [2]), excited by B stars. iHodappI (|l994h observed two 
nebulae in Sh2-233 and noticed a young infrared EC, a s well 
as a probable Herbig-Haro object. IPorras et al.l (|2000|) iden- 
tified two clusters in this region, PCS 2, located around the 
IRAS 05358 + 3543 source and Sh2-233 SE at a separation 
of 1' (0.5 pc) from each ot her for an ado p ted d istance to the 
Sun of 1.8 kpc. However, IChan fc FichI ([l99^ estimated a 
distance of 2.3 ± 0.7 kpc, based on spectral types of stars 
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Figure 6. Same as Fig. |5] for Sh2-235 B, BDSB 73, BDSB 72, 
BDSB 71 and Sh2-232 IR Cluster, respectively. 
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Figure 7. Same as Fig. |4] for the cluster pair PCS 2 and Sh2- 
233SE. 



in the H II region. iPorras et al. I ('2000') detected 92 stars in 
JHK, but the those associated to PCS 2 is 20 and to Sh2-233 
SE is 15 stars. The average extinction was Av = 8.44 ±4.77 
for Sh2-233 SE and Ay = 15.06 ± 3.48 for PCS 2. The es- 
timated ages were 6 Myr for field star s, 3 Myr for S h2-23 3 



SE and less than 2 Myr f or PCS 2. iMao &: Zend l|2004l ). 
based on cluster members (jPorras et al.ll200Gl ). obtained a 
stellar mass of 10.7 Mq for PCS 2 with a star formation effi- 
ciency (SEE) of ~ 4%, and for Sh2-233 SE a mass of 38.4 Mq 
with a SEE of ~ 47%. They suggest that the SEE of PCS 
2 was underestimated because of the upper mass complete- 
ness limit of 1 Mq estimated for PCS 2, and point out that 
the SEE increases with time while the cluster is embedded 
in a molecular cloud. The large SEE difference of the two 
cl usters indicates that Sh2-2 3 3 SE is more evolved than PCS 
2. ISnell. Dickman fc Huan3 (llOQCD identified a CO outflow 
with t wo lobes centred on IRAS 05358 + 3543. Ijiang et al.l 
l|200ll ) named Sh2-233SE and PCS 2 as S233 A and S233 
B, respectively. They argu e that Sh2 - 233SE is less embed- 
ded than PCS 2. Recently. IVan et al.l (|2010l ') derived ages of 
0.3, 0.5, and 1.5 Myr and masses of 45, 30, and 107 M© for 
Sh2-233 SE Cluster, PCS2 and field stars, respectively. 

The nebula G173.58 + 2.45 shows multiple out- 
flow sources, typical of regions presenting youn; 
stellar objects (YSOs) (jWo uterloot fc Brand! 



Shepherd fc Churchwelll 1 19961: [S hcphcrd & Watso^ 



198' 



20021 : 



appears to be associated with a binary system near the 
centre of the G173.58+2.45 Cluster (hereafter G173). We 
discovered an additional EC, CBB 1, at an angular distance 
of 1' from G173. 



3 2MASS PHOTOMETRY 



2MAS^| photometry (Skrutskie et al.ll2006f ) in the J, H and 
Ks bands was extracted in concentric regions centred on 
the coordinates of the ECs (Table [2} using VizieE0. Large 
extraction areas are essential to build RDPs (Sect. [4]) with 
a high resulting contrast relative to the background, and for 
a consistent field star decontamination fSect. [S?T]) . 

3.1 Field-star decontamination 

Eield stars are conspicuous in the observed CMDs of the 
present ECs (Eigs. |4] to [S]!. To uncover the intrinsic CMD 
morphology from the field stars, we apply a field-star decon- 
tamination procedure. The algorithm deals statistically the 
relative number-densities of probable cluster and field stars 



1^ The Two Micron All Sky Survey, 
www. . ipac. caltech. edu/2mass/releases / allsky / 



available at 



Varricatt. Davis fc Adamsonll2005l ). The most extended jet ^ http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=II/246 
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Figure 8. Same as Fig. |4]for tlie cluster pair G173 and CBB 1. 
The diamond is the IRAS 05361+3539 source. 




Figure 9. Colour-colour diagrams for th e decontaminated pho- 
tometry. ISiess. Dufour &: Forestini l|200Cl ) isochrones and redden- 
ing vectors are compared with the PMS distribution. 



in cubic CMD cells that have axes along the J, (J — H) and 
(J — Ks) axes. These are the colours that provide the maxi- 
mum discrimination among CMD sequences for star clusters 
of different ages (e.g. iBonatto. Bica fc Girardll2004 ). 

The algorithm (i) divides the range of magnitude and 
colours of a given CMD into a 3D grid, (ii) computes the 
expected number-density of field stars in each cell based on 
the number of comparison field stars (within la Poisson fluc- 
tuation) with magnitude and colours compatible with those 
of the cell, and (iii) subtracts from each cell a number of 
stars that corresponds to the number-density of field stars 
measured within the same cell in the comparison field. Con- 
sequently, this method is sensitive to local variations in field 
star contamination in magnitude and colours. Cell dimen- 
sions are AJ = 1.0, and A(J ~ H) = A(J - Ks) = 0.2, 
which provide sufficient star-count statistics in individual 
cells and preserve the morphology of the CMD evolution- 
ary sequences. The dimensions of the colour/magnitude cells 
can be subsequently changed so that the total number of 
stars subtracted in the whole cluster area matches the ex- 
pected one {la Poisson fiuctuation). We gave a brief descrip- 
tio n of the field-star decont amin ation procedure. For details 
see iBonatto fc Bical (|2007al ) andU ica. Bonatto fc Camargol 
l|2008l ). 



3.2 EYmdamental parameters 

Usually, decontaminated CMDs of ECs are characterised by 
a poorly-populated MS and a significant number of PMS 
stars. PMS stars are usually located in the same CMD re- 
gion occupied by faint red field stars and consequently one is 
not able to distinguish them, considering only CMD prop- 
erties. In this context, the field-star decontamination pro- 
cedure is essential for the identification of cluster evolu- 
tion sequences. We ado pt isochrones from t he Padova group 
with solar-metallicity (|Girardi et al.l [2002h computed with 
the 2MASS J, H and Ks fihe rs, and the PMS tracks of 
ISiess. Dufour fc Forestiiiil (|2000l ') fitted in J x {J ~ H) and 
Jx {J — Ks) decontaminated CMDs (Figs. |4]to|8} to derive 
the ECs' fundamental parameters. Because of the poorly- 
populated MSs, the 2MASS photometric uncertainties for 
lower sequences, the large population of PMS stars and dif- 
ferential reddening, we decided for the direct comparison 
of isochrones with the decontaminated CMD morphology. 
We carry out eye fits, taking the combined MS and PMS 
stellar distributions as constraint, allowing for differential 
reddening and photometric uncertainties. The cluster sam- 
ple presents a significant fraction of stars redder than the 
youngest PMS isochrone. Probably, the differential redden- 
ing contributed to this (J — Ks) excess towards red colours, 
but part of it can be intrinsic and related to their evolu- 
tionary stage. The youngest stars are the reddest and with 
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Table 2. Derived fundamental parameters for the embedded clusters. 



Cluster 


a(2000) 


(5(2000) 




b 


Av 


Age 


dQ 




Rgc 




(hms) 


(°"') 


n 


n 


(mag) 


Myr 


(kpc) 




(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


KKC 11 


5:41:30 


35:48:49 


173.685 


2.865 


5.5 lb 1.0 


3 lb 2 


2.2 lb 0.5 


9 


.4 ± 0.2 


b bti 784 


5:40:46 


35:55:06 


173.517 


2.794 


4.0 lb 1.0 


3 lb 2 


2.4 lb 0.5 


9 


.6 ± 0.2 


bh2-26o cjZ 


5:41:24 


35:52:21 


1 TO f O A 

173.624 


2.878 


4.0 lb 1.0 


3 lb 2 


2.1 lb 0.5 


9 


.3 ± 0.2 




5:41:07 


DC Af\ or\ 

35:49:30 


T TO r'O A 

173.634 


2.805 


3.8 lb 1.0 


5 lb 2 


2.0 lb 0.6 


9 


.2 ± 0.2 


BDSB 73 


5:40:55 


35:44:08 


173.688 


2.723 


3.8 lb 1.0 


3ib2 


2.1 ±0.5 


9, 


.3 ±0.5 


Sh2-235B CI. 


5:40:53 


35:42:15 


173.712 


2.701 


3.8 lb 1.0 


3ib2 


1.9 lb 0.5 


9, 


.1 ±0.5 


BDSB 72 


5:40:54 


35:40:22 


173.740 


2.687 


3.8 lb 1.0 


3ib2 


2.1 ±0.5 


9, 


.3 ±0.5 


BDSB 71 


5:40:51 


35:38:20 


173.763 


2.660 


3.8 lb 1.0 


3ib2 


2.0 ±0.4 


9, 


.2 ±0.4 


Sh2-232 IR 


5:41:06 


36:09:00 


173.356 


2.973 


5.0 lb 1.2 


3ib2 


1.9 ±0.5 


9, 


.1 ±0.5 


PCS 2 


5:39:13 


35:45:53 


173.481 


2.446 


3.5 lb 1.0 


3ib2 


2.2 ±0.5 


9, 


.4 ±0.5 


Sh2-233 SE 


5:39:10 


35:45:15 


173.484 


2.432 


3.5 lb 0.8 


3ib2 


2.2 ±0.5 


9, 


.4 ±0.5 


G173 CI. 


5:39:28 


35:40:43 


173.581 


2.443 


3.5 lb 0.9 


5ib3 


2.5 ±0.2 


9, 


.8 ±0.2 


CBBl 


5:39:23 


35:41:22 


173.563 


2.434 


3.5 lb 0.9 


5ib3 


2.5 ±0.2 


9, 


.8 ±0.2 


CBB2 


5:41:11 


35:50:10 


173.632 


2.822 


3.8 lb 1.0 


3ib2 


2.0 ±0.6 


9, 


.2 ±0.2 



Table Notes. Cols. 2 to 5: Optimised central coordinates; Col. 6: reddening in the cluster's central region. Col. 7: age, from 2MASS 
photometry. Col. 8: distance to the Sun. Col. 9: Rqc calculated with Rq = 7.2 kpc llBica et al.ll2006l) as the distance of the Sun to the 
Galactic centre. The ECs CBB 1 and CBB2 were discovered in the present work. 



evol ution their colours bec ome bluer and move towards the 
MS (|Lada fc Adamsll 19921 ). Most of these redder stars have 
infrared excesses owing to disks. The isochrone fit gives 
the observed distance modulus (m — M)j and reddening 
E{J — H), which converts to E{B — V) and Ay with the re- 
lations Aj/Av = 0.276, Ah/Av = 0.176, AkJAv = 0.118, 
Aj = 2.76 X E(J - H) and E (J ~ H) = 0.33 x E{B - V) 
l|Dutra. Santiago fc Bicall2002l ). We estimate Av assuming 
a constant total-to-selective absorption ratio Rv = 3.1 and 
adopt the Sun's d istance to the Galactic centre Rq = 7.2 kpc 
l|Bica et al.ll2006l ). The resulting Av, age, and Rqc val- 
ues are given in cols. 4 to 7 of Table [2l respectively. We also 
present in the field subtracted CMDs the reddening vector 
for = to 5. 

In Fig. g] we show the CMDs of the cluster pair Sh2- 
235 {R = 2.5' and R = 1') and CBB 2 (i? = 1'). The aim 
of this procedure is to isolate the population of each EC 
and discard the possibility that we may be working with 
a single object. In the top-panel we present the J x (J — 
Ks) CMDs extracted for Sh2-235 cluster and CBB 2. In the 
middle panels we show the background field corresponding 
to a ring with the same area as the central region. In the 
bottom panels we give the field star decontaminated CMDs. 
For the ensemble we get Av = 3.77±1.0 mag, and a distance 
to the Sun dp = 2.0 d= 0.1 kp c, which agrees with previous 
values (|Georgelin et al.iri973l ). 

Fig. [5] presents the J x (J - Ks) CMDs for FSR 784, 
KKC 11 and Sh2-235E2. These objects are very young clus- 
ter characterised by the presence of a poorly-populated MS 
and a larger number of PMS that appears to be very red- 
dened. These features are absent in the comparison field. 
The gap between the MS and PMS is characteristic of young 



clusters. The combination of MS and PMS shows that these 
objects are ECs. 

We show in Fig. [6] the CMDs of the ECs Sh2-235 B, 
BDSB 73, BDSB 72, BDSB 71 and Sh2-232 IR Cluster. 
These objects present a poorly-populated MS and a rather 
populous PMS. 

Figs. [3 and [g] show the CMDs of cluster pairs. Fig. [3 
shows the CMDs of the pair PCS 2 (7? = 1') and Sh2-233 
SE Cluster {R = 2' and R = 1'). Fig. [8]presents the CMDs 
of the pair CBB 1 (R = 1') and G173.58 Cluster (R = 2' 
and R — 1'). Both objects present CMDs characteristic of 
very young clusters. The pair separations are small and have 
a number of stars in common. We do not discard the possi- 
bility that these two objects will merge. The — 2.2 ± 0.5 
kpc for Sh2-233 SE agrees with d© — 2.3 ± 0.7 estimated by 
IChan fc FichI ^9^)- The diamond in the CMD of the G173 
is the IRAS 05361+3539 source (Sect. [2]), which appears to 
be a PMS star approaching the MS. 

We estimate Av for BDSB 71 and Sh2-233SE Clus- 
ter assuming a constant total-to-selective absorption ra- 
tio Rv = 5.0, which is t ypical of dense clouds 
(|Cardelh. Clayton fc Mathislll989l '). BDSB 71 presents Av = 
6.0 ± 1.6 mag for a distance from the Sun of d© = 1.3 ± 0.7 
kpc. For Sh2-233SE Cluster the values are Av = 6.4 ± 1.6 
mag and d© = 1.6 ± 0.5 kpc. 

Fig. llOl shows the number-density distribution for prob- 
able stars in the decontaminated J x ( J — Ks) CMDs. The 
densest parts of those Hess diagrams correspond to the PMS. 

The present sample of ECs suggests some age spread 
implied mainly by the PMS stars. This suggests some se- 
quential star formation (Fig. Ilip . However, the visual ex- 
tinction of the same ECs are much smaller than the values 
observed for individual stars in these clusters which suggests 
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Figure 10. Number-density (stars arcmin ^) distribution for probable stars present in the decontaminated CMDs of FSR 784 (top 
left), Sh2-235 Cluster (top right), Sh2-235E2 (bottom left) and Sh2-235B (bottom right) respectively. 



that these stars are probably embedded in dense structures, 
which explains the larger values of Av estimated in the pre- 
vious works, especially in the case of individual star analyses. 



Figs. [2l l3l and [TTI together with the age derived for the 
Sh2-235 Cluster (« 5 Myr), indicate that the O star (si 1 
Myr) is not a member of this cluster. A possible explana- 
tion for the presence of the O star near this cluster is that 
winds from Sh2-235 Cluster, colliding with the surrounding 
gas, might have originated this star, CBB 2 and other clus- 
ters in the neighborhood. Assuming v ~ 20 km for th e 
dense gas involved in the expansion (|Kirsanova at alj|2008l ). 
the Sh2-235 Cluster might be responsible for sequential star 
formation across a region of radius ~ 10 pc (Fig. Ilip . which 
includes also the clusters in a row in the southwest direc- 
tion. Sequential star formation is possible also for the pairs 
G173 and CBB 1, and Sh2-233SE Clu ster and PCS 2. Re- 
cently, iDewangan fc Anandaraol (|201ll ) point out that star 
formation continues to occur in the Sh2-235 complex, mainly 
within the ECs. They identified 86 Class O/I and 144 Class 
II YSOs, which reinforce the possibility of a sequential star 
formation event. 



3.3 Colour-colour diagrams 

Colour-colour diagrams are useful tools to investigate the 
nature of ECs. We show in Fig. [5]the decontaminated near- 
IR colour-colour diagram (J — Ka) x (H — Kb) of the member 
stars , together with PMS tracks l|Siess. Dufour &: Forestinil 
l2000l ). set with the reddening values derived above, to es- 
timate ages. As a consequence of the presence of the PMS 
stars in the cluster, it is expected that some stars present 
near-IR excess. As expected from the CMDs of ECs (Figs.|3] 
to [8]), a significant fraction of the stars appears to be very 
reddened. Most stars, specially MS stars, have [H — Kg) 
colours close to the isochrone, within the uncertainties. Be- 
sides, most of the very red PMS stars are displaced parallel 
to the respective reddening vectors. However, a significant 
fraction appears to present an abnormal excess in (J — Ks) 
and {H — Ks), especially Sh2-235B, which may come from 
PMS stars still bearing circumstellar discs. MS stars lie on 
the blue side of the diagrams and there occurs a gap between 
MS and PMS stars in the CMDs. 



4 CLUSTER STRUCTURE 

The structure of the ECs is analysed by means of the stellar 
radial density profile (RDP), defined as the projected num- 
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Figure 11. Schematic distribution of cluster positions and ages. 
Asterisks are O stars, and the plus sign is a B star. It also shows 
the nebulae Sh2-235 and Sh2-231 as reference. 



ber of stars per area surrounding the cluster centre. RDPs 
are built with stars selected after applying the respective 
colour magnitude (CM) filter to the observed photometry. 
CM filters isolate the probable cluster sequences excluding 
stars with colours different from those of the cluster se- 
quences (e.g. , Bonatto fc Bica 2 007ai . and references therein). 
However, residual field stars with colours similar to those 
of the cluster are expected to remain inside the CM filter. 
They affect the intrinsic stellar radial distribution profile in 
a degree that depends on the relative densities of field and 
cluster. The contribution of these residual field stars to the 
RDPs is statistically quantified by means of comparison to 
the field. In practical terms, the use of the CM filters in clus- 
ter sequences enhances the contrast of the RDP with respect 
to the stellar field. The CM filters are shown in Figs. |4]to 
[8] as the shaded area superimposed on the field-star decon- 
taminated CMDs. 

To minimise oversampling near the centre and under- 
sampling for large radii, the RDPs are built by counting 
stars in concentric rings of increasing width with distance 
to the centre. The selected number and width of rings pro- 
duce RDPs with adequate spatial resolution and moderate 
1(7 Poisson errors. The residual background level of each 
RDP corresponds to the average number of CM-filtered stars 
measured in the comparison field. 

The cluster structure was derived by mea ns of a King- 
like profile, which is similar to a two-parameter [Kind (|l962 ) 
model that describes the intermediate and central regions of 
globular clusters. The fit was performed with a non-linear 
least-squares routine that uses the errors as weights. The 
best-fit solutions are shown in Fig. [12] for 3 relatively popu- 
lous ECs as a solid line superimposed on the RDPs. The pro- 



file is expressed as (t{R) = (Tbg +croK /{! + (R/ Rcore) , where 
ffbg is the stellar background surface density, aoK is the cen- 
tral density relative to the background level and Rcore is 
the core radius. The cluster radius (Rrdp) and uncertainty 
can be estimated by considering the RDP fluctuations with 
respect to the residual field. Rrdp is the distance from the 
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Figure 12. Stellar radial densities profile (filled circles) built witfi 
colour-magnitude filtered photometry, centred in the coordinates 
of Sh2-235 Cluster, FSR 784 and Sh2-235E2. Solid line: best-fit 
King profile. Horizontal shaded region: stellar background level 
measured in the comparison field. Gray regions: Icr King fit un- 
certainty. 



cluster centre where RDP and comparison field become sta- 
tistically indistinguishable. Small variations in the RDPs are 
probably due to the presence of other clusters and/or en- 
hanced dust absorption. The derived structural parameters 
are given in Table [3] 

In Fig. [12] we show the RDPs of Sh2-235 Cluster, Sh2- 
235 East2, and FSR 784. For Sh2-235 Cluster and Sh2-235 
East2. Overdensities show up in their RDPs. 

Fig. [13] shows the RDPs of Sh2-235B Cluster, BDSB 
71, BDSB 72, BDSB 73, Sh2-232 IR Cluster, KKCU, CBB 
2 and the pairs CBB 1 and G 173, PCS 2 and Sh2-233 SE. 
The RDPs are typical of ECs of low mass and/or initial 
evol utionary phases, a nd cannot be fitted by King's pro- 
file (|Soares et al.l |2005| ) . They present bumps and dips as 
compared to field stars. The depression in star counts in 
the central region of some ECs is possibly due to strong 
dust absorption, crowding or structured cores. Sh2-232 IR 
Cluster has a high cluster/background density contrast, but 
the profile is irregular. The Sh2-235B cluster profile includes 
several of the small clusters nearby like BDSB 71, 72 and 
73. A deeper photometry is required for the analysis of the 
structure of these objects. 

The structure of young populous ECs can be gener- 
ally characterised by an RDP with multiple peaks on a 
large spatial scale or centrally condensed with an RDP that 
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Table 3. Structural parameters. 



Cluster 


(1') 


O"0A' 


Rcore 


Rrdp 




Rcore 


Rrdp 




Ai? 




(pc) 


(* pc^^) 


{pc) 


(pc) 




(') 


{') 




(') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


FSR 784 


0.69 


56.7 ±4.0 


0.25 ±0.01 


2.1 ±0.7 


27.02 ± 1.9 


0.36 ± 0.02 


3.0 ±1.0 


20 


- 60 


Sh2-235 E2 


0.60 


71.7 ±25.0 


0.13 ±0.03 


1.2 ±0.3 


25.8 ±9.0 


0.21 ± 0.05 


2.0 ±0.5 


20 


- 80 


Sh2-235 CI. 


0.56 


138.1 ±27.4 


0.10 ±0.01 


1.4 ±0.3 


43.32 ±8.6 


0.18 ± 0.02 


2.5 ±0.5 


20 


- 80 



Table Notes. Col. 2: arcmin to parsec scale. To minimise degrees of freedom in RDP fits with the King-like profile (see text), tTi,g was 
kept fixed (measured in the respective comparison fields) while (Tq and Rcore were allowed to vary. Col. 9: comparison field ring. 



Table 4. Stellar mass estimate for more populous ECs 



MS PMS MS + PMS 



Cluster AmMS N M N M N M 

(Mq) (stars) (Mq) (stars) (Mq) (stars) (Mq) 

(1) (2) (3) (4) (5) (6) (7) (8) 

FSR 784 2.90-4.30 3±1 9 ±4 65 ± 15 39 ± 9 68 ± 16 48 ± 13 

Sh2-235E2 2.50-5.25 6±2 24 ± 9 41 ± 6 25 ± 4 47 ± 8 49 ± 13 

Sh2-235C1. 1.30-6.25 5±2 22 ± 7 67 ± 15 40 ± 9 72 ± 17 62 ± 16 



Table Notes. Col. 2: MS mass range. Cols. 3-6: stellar content of the MS and PMS stars. Cols. 7-8; total (MS+PMS) stellar content. 



can be described by a King's law (|Lada fc Ladal [20031 ) . al- 
though the present objects are not virialised. However, the 
differential dust absorption produces conspicuous variations 
in the RDPs of the present objects. King profile describes 
the structure of clusters close to spherical symmetry and 
centrally concentrated. However, many young clusters are 
substructured or asymmetric, deviating significantly from 
this shape, and therefore cannot be fitted by King's law 
ICartwright fc Whitwo"rthll2004 1 Gutermuth et al. I [2005h . In 
Fig. 1141 we show the spatial distribution of stars in the de- 
contaminated photometry of the three ECs that follow a 
King-like profile (FSR 784, Sh2-235 Cluster and Sh2-235E2) 
and two representative cases of objects that do not (KKC 11 
and BDSB 73). The former are centrally concentrated and 
nearly circularly symmetric. The cavities and overdensities 
in the stellar distribution (Fig. I14|) can be seen as bumps 
and dips in the RDP (Fig. 1111) . On the other hand, objects 
like BDSB 73 that are not centrally condensed and KKC 11 
with more elongated shape do not follow a King profile. The 
multiple peaks in the RDPs of these ECs may be a fractal 
effect. If these objects survive the primordial gas expulsion, 
they may undergo merging evolving into a relatively smooth 
structure. 

The angular distribution of decontaminated stars 
(Fig. ll4p . used in the CMD construction, reproduces the dis- 
tribution of stars in the RDPs built with filtered photometry 
(Figs. [T2I and [T3)l . supporting the consistency of our results. 



5 MASS ESTIMATES 

Given the poorly-populated nature of the MS, we simply 
counted stars in the CMDs (within the region R < Rrdp), 



and summed their masses as estimated from the mass- 
luminosity relation implied by the respective isochrone so- 
lution (Sect. |3^ . The results are given in Table |4l 

All the ECs clearly present distinct populations of MS 
and PMS stars (Figs.|4]to|8]). However, given the differential 
reddening, it is not possible to attribute a precise mass value 
for each PMS star. Thus, we simply count the number of 
PMS stars and adopt an average mass value for the PMS 
stars to estimate upms and mpMS- Assuming that the mass 
distribution of the PMS stars also follows Kroupa's (2001) 
MF, the zsc, average PMS mass - for masses within the 
range 0.08 < m(MQ) < 7 - is < tupms >« 0.6Mq. Thus, 
we simply multiply the number of PMS stars (Table 8) by 
this value to estimate the PMS mass. Finally, we add the 
latter value to the MS mass to obtain an estimate of the 
total stellar mass. These values should be taken as lower 
limits. 



6 RELATIONS AMONG ASTROPHYSICAL 
PARAMETERS 

N-body simulations of massive star clusters that include the 
effect of gas removal (e.g. Goodwin & Bastian 2006) show 
that the phase of dramatic core radii increase may last about 
10-30 Myr. Mass segregation may also lead to a phase of core 
contraction, with high mass stars more concentrated in the 
core while low mass stars are transferred to outer parts of the 
cluster. In this context, we suggest that most objects of our 
sample have not yet removed completely their primordial 
gas and are not in this expansion phase. As a consequence 
the core and cluster radii remain small. However, these ECs 
may be intrinsically small as well. 
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Table 5. Integrated colours and magnitudes. 



Magnitude 



Colour 







Apparent 






Absolute 




Reddening 


Corrected 


Cluster 


J 


H 




J 


H 




(J~H) 


(J-K,) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


FSR 784 


12.9 


9.7 


8.8 


-0.02 


-2.8 


-4.0 


2.82 ±0.14 


3.57 ± 0.14 


Sh2-235 E2 


13.2 


11.7 


10.6 


0.5 


-0.6 


-1.5 


1.11 ± 0.33 


1.99 ± 0.31 


Sh2-235 CI. 


8.5 


8.2 


7.9 


-4.0 


-4.0 


-4.0 


-0.05 ±0.04 


0.01 ± 0.04 



Table Notes. Col. 2-4: apparent magnitude. Cols. 5-7: absolute magnitude. Cols. 8-9: (J — H) and (J — Kg) colours. 
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Figure 13. Stellar RDPs for the remaining ECs built with colour- 
magnitude filtered photometry. 



iLada fc Ladal l|2003D argue that the structure in ECs 
reflects the underlying structure in the dense molecular gas 
from which they formed and suggest two possible structures 
for these objects: those similar to classical OCs (radially 
concentrated that fit a King law) and those that exhibit a 
density profile with multiple peaks. The latter kind of ECs 
present a fractal-like structure with smaller substructures 
or mini clusters that are probably a consequence of the 
fact that these o bjects are formed in a CMC with a frac- 
tal substructure (ISchmeia et al ] l2008l : iLomax et all l2010l : 
I Sanchez et all |2010| . and references therein). They evolve 
with time and their fates depend on the processes of dis- 
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Figure 14. Angular distribution of the decontaminated stellar 
content for FSR 784, Sh2-235, Sh2-235E2, KKC 11 and BDSB 
73. BD+35°1201 is indicated as an open circle around the star. 



solution, which they will undergo. If the ECs do not endure 
the action of the dissolution processes, they might evolve 
into a homogeneous distribution of stars and eventually dis- 
perse. Otherwise their fate will be a centrally concentrated 
distribution of stars or a bound cluster. This is possible for 
the lower mass ECs and cluster pairs or multiples (Figs.|3]to 
[H|. It would also be important to have deeper observations 
of the less-populated clusters with irregular RDPs (Fig. llSf) . 
Probing fainter PMS stars might contribute to the construc- 
tion of better sampled stellar density profiles, like those of 
the intrinsically populous and massive clusters in the present 
star-forming complex (Fig. [HJ. 

When the projected mass density of a star cluster fol- 
lows a King- like profile (iBonatto fc Bical l2008ah . the clus- 
ter mass (Mcius) can be expressed as a function of the 
core radius and the central surface mass-density (a mo) ac - 
cording to M^ius « IS.Sctmo Rc (|Bonatto fc Bical 120091 '). 
Fig. [15] (top-right panel) shows the distribution of our pop- 
ulous ECs in the plane core radius (Sect. \Q vs cluster 
mass (Sect. [Sjl. Clearly, our ECs (together with the refer- 
ence sample) distribute parallel to the above relation, be- 
ing constrained within King-like distributions with 15 < 
o"mo(Mqpc~^) < 600 (these limits take the uncertain- 
ties into account). This correlation between cluster mass 
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Figure 15. Present cluster properties compared to literature 
ones. Top-left panel: a ge and Rcore r elatio ns. The dashed line 
is the logarithmic fit of iBastian et al] | |200S| ) {Rcore{pc) = 0.6 X 
ln{age[Myr]) — 0.25) for M51 clusters. The solid line is our fit for 
Galactic clusters. Bottom- left panel: IzgcI for clusters of the top 
panel. Top right: core radius and cluster mass follow the relation 
Mcius = 13.8(TM0^ci with varying values of ctmo'i reference sam- 
ple OCs are shown as gray circles. Bottom right: cluster density 
vs. radius. Arrows indicate lower limits to the cluster mass and 
density. 



7 CONCLUDING REMARKS 

In the present work we performed a field-star decontami- 
nated 2MASS analysis of 14 ECs in the H II regions Sh2-235, 
Sh2-233, Sh2-232 and Sh2-231, and other small nebulae in 
the area. We were able to derive astrophysical parameters 
and investigate the nature of young clusters embedded in 
these H II regions. Fundamental parameters were obtained 
for all ECs analysed, but structural parameters were de- 
rived for FSR 784, Sh2-235E2 and Sh2-235 Cluster. The de- 
contaminated CMDs exhibit a poorly-populated MS and a 
large fraction of PMS stars, and suggest some age spread, 
as expected for sequential star formation process. Two new 
ECs (CBB 1 and CBB 2) were discovered in this work. 

The age, size and location of KKC 11, FSR 784, Sh2-235 
Cluster, CBB 2, and Sh2-235E2 are consistent with the col- 
lect and collapse scenario. The CMDs of Sh2-235 B, BDSB 
71, BDSB 72, BDSB 73, Sh2-232 IR Cluster, G173, CBB 
1, and PCS 2 suggest that these objects are ECs. The en- 
hanced dust absorption mainly in the innermost region in- 
dicates that deep photometry is necessary to derive their 
structural parameters. 

The present ECs have core and cluster radii smaller 
than clusters at the same Galactocentric distance and age. 
Probably most of them have not expelled the primordial gas 
completely, and thus have not expanded. 
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and core ra dius is consistent with the mass-radius rela tion 
suggested bv lPortegies Zwart. McMillan fc Gielej (|20ld ) for 
massive cluste rs younger than 100 Myr and ext ended for less 
massive ones (|Camargo. Bonatto fc BicalbOlOl ). 

Since, the cluster radius and mass were estimated, 
we compute the cluster mass density p{Mq pc~"^) = 
^McIusRrdp- The results are shown in the plane 
Rrdp vs. p (Fig. 1151 bottom-right panel). Despite the er- 
ror bars, the density decreases with cluster radius as p oc 
-Rpnp^°^\ similarly to t he cl uster sample studied by 
ICamargo. Bonatto fc Bical l|2010l ). We also show in Fig. [12] 
the relation between age and both Rcore and \Zgc\ for the 
present ECs. We fit an empirical curve Rdjpc) = 0.27 x 
In [age[Myr\) — (the solid line) to Galactic OCs younger 
than 100 Myr . The dashed line is the observed relation for 
M51 clusters (|Bastian et al.|[2008t ). The present clusters be- 
have as relatively low mass ECs, as expected. 
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